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Abstract 


In this report, a new method of thermal protection for transatmospheric vehicles is 
introduced. The method involves the combination of radiation, ablation and transpiration 
cooling. By placing an ablating material behind a fixed-shape, porous outer shield, the 
effectiveness of transpiration cooling is made possible while retaining the simplicity of a 
passive mechanism. A simplified one- dimensional approach is used to derive the governing 
equations. Reduction of these equations to non-dimensional form yields two parameters 
which characterize the thermal protection effectiveness of the shield and ablator combina- 
tion for a given trajectory. 

The non-dimensional equations axe solved numerically for a sample trajectory corre- 
sponding to glide re-entry. Four typical ablators are tested and compared with results 
obtained by using the thermal properties of water. For the present level of analysis, the 
numerical computations adequately support the analytical model. 
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Nomenclature 


A 
c 
F 
G 
Ho 
HC 
HL, L 
m 
Qo 
Q r 

Q 

R 

t 

T 

V 

x,y,z 


dimensionless ablation parameter defined in text 
heat capacity (specific heat, constant pressure) J/kg/K 
dimensionless heating rate function 
dimensionless total heat input function 
total (stagnation point) enthalpy of free stream, J /kg 
threshold enthalpy of ablating material, J /kg 
enthalpy of ablation, J /kg 
mass per unit area, kg/m 2 

stagnation point convective heating rate, W/m 2 
re-radiation heating rate, W/m 2 

total heat input (heating rate integrated over time), J/m 2 
dimensionless radiation parameter 
time, s 

temperature, K 
dimensionless mass flux 
cartesian coordinate system 


Greek Symbols 


a absorptivity coefficient 

/? parameter, used in reduced form of equations in text 
6 dimensionless integral thermal thickness 
e emissivity 

rj transpiration coefficient 
p density (kg/m 2 ) 

a Stefan-Bolzman constant = 5.67 x 10 ~ 3 W/m 2 K 4 


n 


k conductivity, W/m/K 


Subscripts 

a ablator property at time of ablation 
b solid ablation material, virgin state 
o stagnation point value 
r radiative 

s ablation surface prior to ablation 
w shield wall 
oo free-stream property 

Superscripts 

* critical, maximum value 
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1.0 Introduction 


Radiating ablation materials have seen extensive use for the protection of aerospace 
vehicles over the past thirty years, notable examples being the ICBM nose cones and the 
manned spacecraft Mercury, Gemini and Apollo. These materials have successfully served 
the purpose of thermal protection but have distinct disadvantages for reusable vehicles, 
particularly those for which the gometry determines the aerodynamic performance. The 
use of such ablation materials in regions of localized high heating rates, such as the nose 
or leading edge of the wing leads to significant changes in shape over critical aerodynamic 
surfaces with consequent adverse effects on performance characteristics. 

In order to overcome such disadvantages while retaining the inherent passive simplicity 
of the radiating ablation approach, a new concept is studied in this paper: an internal 
ablator protected by a fixed-shape porous outer radiating structure. This approach permits 
leading edges and other regions of intense localized heating to be radiatively cooled over 
most of the flight trajectory. When excessive heating occurs however, ablation results 
and causes the transpiration of gaseous products through the outer radiating structure, 
thus controlling its temperature. The introduction of ablation products into the external 
boundary layer reduces the net heat flux to the outer surface thus providing a passive 
thermal protection mechanism. 

This approach to thermal protection depends on the judicious choice of thermal prop- 
erties of the outer structure and the internal ablator. The analysis given in this paper 
provides the thermal parameters which assure that the external surface temperature of the 
radiating structure will not exceed specified limits during the heating period while avoiding 
excessive mass loss from the internal ablation material. Approximate methods of analysis 
are used to identify the important thermal parameters in dimensionless form which may 
be used for the selection of material properties. 
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1.1 Description of the Concept 


The radiation/ablation concept considered here consists of a porous fixed-shape outer 
structure, the temperature of which is limited to a value T* (see Fig. 1.1). Behind 
this outer structure is placed an ablator of sufficient thickness such that its back surface 
temperature T\, remains constant while its front surface when heated will ablate at a 
characteristic temperature T a . During ablation gaseous products transpire through the 
fixed outer structure into the external boundary layer thereby providing the means for 
reducing the heat transfer rate to the outer surface. Heat transfer between the outer 
structure and the ablation surface is by radiation across the gap, thermal conduction in the 
gaseous products within the gap is ignored. A planar geometry is considered corresponding 
to the region near the stagnation point of the vehicle nose or near the stagnation line of 
the wing. Extension of the concept to other geometries should be straightforward. 

2.0 Concept Analysis 

In the present analysis it is assumed that the outer structure is of high thermal con- 
ductivity and attains a uniform temperature T w {t) which must be limited throughout the 
flight to a value less than some maximum value T*. This structure is perforated to allow 
for the transpiration of ablation products into the boundary layer which thereby provide 
a reduction in the heating rate. The outer structure radiates outward at a rate q r (t) and 
also radiates inward at a rate aq r (t ) to the ablation surface, where a is the absorption 
coefficient of the ablation material. 

It is assumed that the ablation material has a surface temperature T s (t) which increases 
to a fixed value T 0 when ablation occurs. The ablation material is considered sufficiently 
thick to permit its back surface temperature to maintain a constant value Tj. 
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Thermal Balance Equations 


In the absence of ablation the thermal balance at the outer surface can be written 


(1 + a)q,(t) = ,„(() (2.1) 

where q T (t ) = creT* and q 0 (t) is the convective heating rate for the given trajectory. When 
ablation occurs, however, transpiration through the outer structure into the boundary layer 
flow causes a reduction in the heating rate, which can be expressed in simplified form as 

(1 + £»),,(<) = - r/H^ (2.2) 

In equation (2.2) H a is the total stream enthalpy, r/ is an efficiency factor which depends 
on the state of the boundary layer (laminar or turbulent) and on the molecular weight of 
the ablation products, and dm/dt is the rate of mass transpiration. The factor r/ lies in 
the range 1/3 < rj < 2/3 for typical ablation materials. The quantities q r (t ) and dm/dt 
are in general unknown and must be determined by considering the heat radiated to the 
ablation surface and the conduction of heat within the ablation material. 

It is convenient to define a thermal integral thickness 0(t ) within the ablation material 
such that the heat content of the material is pc(T s — T b )9 where (T s — T b ) is the temperature 
rise at the surfacef 1-3 ^ Thus 


= /‘ 


0 T -T b 


6 = I ^ — -±dz 

-00 T 3 — T& 


(2.3) 


The overall heat balance, in differential form, for the ablation material, is then 


cq r (t) = [c,(T 0 - 7J) + L) ~ + *4 |(T, - T t )6] (2.4) 
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where L is the enthalpy of ablation (vaporization) and c b {T a — T b ) is the threshold enthalpy 
rise prior to vaporization. An additional equation representing the heating rate boundary 
condition at the ablation surface is written 


“**> =L % + h {£), M <2 ' 5) 
Equation (2.5) represents a heat transfer rate balance. Since the ablator is considered to 
have thickness much larger than 6 the approximation 

, fdT\ k b (T,-T b ) 

h {n) 0 = —e— 

can be made, with the result that equation (2.5) is written 


a ir (t) = L— + 


(2.5a) 


Finally the unknown quantity q T (t ) can be eliminated from equations (2.4) and (2.5a) 
through the use of equation (2.2) to give 


(iTa) q °^ = h (T “ - T >) + 1 + (iqrj) i H ° 


^ + Pi [(T, - TM (2.6) 


and 




1 dm k b (T s - T b ) 


d 


(2.7) 


where q 0 (t) is the known heating rate that would occur at the outer surface in the absence 
of ablation and H 0 {t) is the total stream enthalpy, both determined from the vehicle 
trajectory. 
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Equations (2.6) and (2.7) must be solved either for T s and 9 during periods of non- 
ablation when dmjdt = 0 and for dm/dt and 9 during periods of ablation, when T a = T a . 

2.1 Pre-Ablation Heating 

During the pre-ablation heating period, dm/dt = 0 and equation (2.6) may be inte- 
grated to give 


(iT^) = pbCb{T> ~ Tb) 6 

while equation (2.7) reduces to 


( 2 . 8 ) 


( oc \ h(T,-T t ) 

liTsJ q ° (t) ? 

where the total heat input is 

<?.«= /'«.M 

Jo 

Equations (2.8) and (2.9) give the solution 


and 


( a ) 


\l+aj 

pbCbh 


( 2 . 10 ) 


9 = 


Qo(t)h 


[q 0 (t)pbCb \ 


( 2 . 11 ) 


giving the temperature at the ablator surface and the heat penetration thickness 6. From 
equation (2.10) the time t a at which the surface temperature reaches the value T a is given 


5 


implicity by 


Qo(^o)?o(^o) — QaQ a — Pb c b^b(^a — TJ,) 2 (2.12) 

and the thermal thickness at this time is O a = [Q a ^bl %Pb^b[ z - Prior to ablation, the surface 
temperature of the outer structure T w is given by equation (2.1) 


aeT^ = q T (t) = 


go(0 

1 + a 


or 


T w (t) = 


go (t) 


(2.13) 


L(l + a)creJ 

If the maximum outer surface temperature permitted is T* (corresponding to a heating 
rate q* = aeT* 4 ) then the time t a when ablation occurs must satisfy the implicit condition 


go(^a) < 1 
(1 + a)q; ~ 


(2.14) 


This must be achieved through proper choice of the ablation material properties 
a, /?(,, kb, and T a as seen from equation (2.12) 


2.2 Ablation Heating 

During the period of ablation, for which T s = T a , equations (2.6) and (2.7) become 


(its) - h T - - 1 + fe) 3? + - r ‘)f (*•«) 
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and 




1 dm k b (T a - T b ) 
dt + 9 


(2.16) 


which may be solved for dm/dt and 0. These equations are solved in dimensionless form 
in Appendix A and approximate analytic solutions are given under the assumption that 
the enthalpy ratio: 


Cp(T a - T h ) 

C P (T. - Tb) + (*) qH. 


(2.17) 


In view of the large value of H 0 compared with all other terms in (2.17) this assumption 
is generally valid. 

The approximate solution is written: 


d/6 a = [1 + 2(3 a (t/t a - 1)]* (2.18) 

and from 2.15 


dm fa} 

«o(*) 

1 

+ 

rH 

1 

c b{T a — Tfc) + L + 

(ife) 

q 0 [l + 2/3 a (t/t a -l)]$. 


(2.19) 


where fl a = q a t a /Q a . The approximate solutions given by (2.18) and (2.19) are valid only 
as long as the inequality (2.17) holds. 

During ablation the temperature of the outer structure is found from 

fZkV _«L 

\T*J q • 
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where q T is given by equation (2.5a) with 6 and dm/dt determined from equations (2.15) 
and (2.16) respectively. These simplified analytical solutions jure presented for the purpose 
of deriving the relation (2.21) below. The result is written 


( ZA Qa . f Qo Qg \ ^ 

It*; “ (l + a)q; [d U e )' L + (^) qHo 


( 2 . 20 ) 


The solution 8/d a is derived in Appendix A, and the approximate solution is given 
by equation (2.18). In order to assure that T w /T* < 1 during ablation the latent heat 
of ablation, L, must be such as to allow the release of an adequate amount of ablation 
products to cool the outer structure. Equation (2.20) gives T w /T * < 1 if 


?a 


(1 + a)q: 
or, after rearrangement of terms: 


_ , (s* _ M L 

;[e W */> + (* 


J 


< 1 


I ^ f [go/gq ~ C 1 + Q )gr/ga] 

L > \[(1 + a)q;/q a -e a /6}.[(r£;)T ] H 0 } 

Thus, the ablation material must satisfy the inequality (2.21) to maintain the shield tem- 
perature at or below the required maximum. 

2.3 Post-Ablation Heating 

From equation (2.16) it can be seen that ablation ceases ( dm/dt = 0) when 


1 . 


( 2 . 21 ) 



h(T. - T t ) 

«((,) 


( 2 . 22 ) 


So that t = t e satisfies the relation q e /qa = 0 o /0e where q e — q 0 (t e ) and 6 e = 8(t e ). 
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For t = t e , equations (2.6) and (2.7) with dm/dt = 0 are again used to determine 9 
and T, — T b . Equation (2.6) is integrated to give 

(iTs) qM)ii = P>C>(T - ~ Tt)e ~ " Ti)e ' (2 - 23) 

and 



h(T, - Tj,) 
0 


where q 0 (t),t e ,q e ,T a and T b are all known. Eliminating 9 gives 


(2.24) 


Pb c b k h {Ts - T b f = [(~L^) (Q 0 - Q e ) + p b c b (T a - T b )9 e ] (---) q Q (2.25) 

where the change in total heat input is given by integral 

/%o(CK = Qo(t) - Qo(te) = Qo-Q e 

Jte 

or alternatively, by substituting for 9 e from (2.22) 


pmh(T, - Tbf = (^) ’(<?„- Q,)q, + p*MT. - Tbf (2.25a) 

which may be used to solve for the ablation surface temperature T, during the post- 
ablation period. Similarly, elimimating (T, — T b ) from equations (2.23) and (2.24), the 
thermal thickness 8 for i > t e is given by 


8{t > t e ) = 


’ (Qo-Qe) 1 

Qo 



PbChhjTg — T;,) 2 

?o?e 



i. 


2 


(2.26) 
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and the surface temperature of the outer structure is again given by equation (2.13) i.e., 


or, in dimensionless form, 


T w = 


Qo(t) 


(1 + a)<xe 


(2.27) 


T* 


?o(0 


1 * 


(2.27a) 


U 1 + a )9r*J 

Thus, the history of outer structure temperature, T w , the ablation surface temperature, T s , 
the ablation rate, and thermal penetration thickness into the ablator may be determined 
in terms of the material thermal properties and the imposed heating rate q 0 (t). 

3.0 Non-Dimensional Analysis 

The numerical solution of the governing equations introduced in the previous section 
is simplified by reducing them to non-dimensional form: Hence, it is convenient to char- 
acterize the heating rate history q 0 {t) in the form 


9o(0 = q m F(t/t m ) (3.1) 

where q m is the maximum heating rate during the flight trajectory, occurring at time t m . 
Similarity the total heat input Q 0 (t) is expressed as 


Qo(t) = QmG(t/t m ) 


(3.2) 


The dimensionless total heat input function G is given by the ratio of integrals 


C(i/i 1 — x/ m F(tjtrn)d{t/i m) 

{/m) - f, 1 F(t/t m )d(t/t m ) 


(3.3) 
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Thus, 0 < F < 1 and 0 < G < 1 in the interval t <t m . In the pre-ablation period 
the temperature at the ablation surface, given by equation (2.10) becomes 

( T,-T b \ (FG)i 
\T a -T b J A? 

where the ablation parameter introduced is defined by 

1 + Q \ 2 PbCbkbi^a — %) 2 
® ' QmQm 

The ablation parameter is a dimensionless number associated with the ability of the ablator 
to conduct heat when subject to the imposed heating rate q m . The thermal penetration 
thickness, equation (2.11), is written 



(3.4) 


6 = 



(3.5) 


and at t = t a 


e a = Qmh 1 * 

The surface temperature of the outer structure is now given by 


(3.6) 


T 

•L W 

T* 




li 


U 1 + <*)£ J 


[F(t/tm) ]* 


Thus in order to assure that T w < T* we must have the inequality 


(3.7) 


F(tjtm ) < 


(1 + a)g; 

9m 


= R 


(3.8) 
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The radiation parameter introduced, R = (1 + a)q*/q m , is a dimensionless number 
associated with the ability of the outer surface to reject heat by radiation. The heating 
rate (and, indirectly, the time t a ) at which the ablation surface reaches the temperature T a 
is also a function only of the ablation parameter A and satisfies the condition, from (2.12), 


F(tJt m )G(t a /t m ) = A 


If G is first expressed as a function of A, this equation may be solved to give 


(3.9) 


F(t a /t m ) = f(A) 

Thus in order to assure that Tw < T*, we must have 

f(A) < R (3.10) 

This important relation delineates a characteristic curve relating the ablation and 
radiation parameters. For a given trajectory, the range of values for these parameters is 
separated into two regions, one for which the shield/ablator combination is an acceptable 
design and another for which the combination will not protect the shield material from 
excessive heating. The shape of the curve given by (3.10) will generally vary for different 
trajectories, but the essential relation will remain, namely that the curve will separate the 
range of values for the ablation and radiation parameters into the two regions described 
above. This will be discussed in greater detail in the next section. 

During the ablation period the solution depends on the parameter /?„ = q a t a /Qa which 
is now rewritten as: 




( Qm tm \ 

'F(t a /t m )t a /t m 

\Qm) 

G(t a /t m ) 
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where F(t a /t m )t a /t m /G(t a /t m ) = F a t a /t m /G a is also a function of A. During ablation the 
outer surface temperature is given by equation (2.20) which may be written 


Iz. = lEi\ e -±+(L- e -±\ L l ! 7 

T * 9) ' L + ariH 0 /(l + a)\j 

The condition on the heat of ablation L, equation (2.21), then becomes 


(3.11) 


1 _ 

L 


> 


{ 

1 

cs 

1 

IV 

1 

i 

1 

11 

[R-(6J9)F a \ 

ar]H 0 /(l + a) 


max 


(3.12) 


This condition assures that ablation will occur as necessary to avoid excessive temperatures 
at the outer surface. 

The governing equations during the ablation period may be simplified by making the 
integral thermal thickness and mass flux non-dimensional. The dimensionless quantities 
are: 


^(Thermal thickness) = 9/9 a 


and 


V (mass flux) = 


fo(T. Z £6 )+£+ (gfefo H 0 ]dm/dt 
aq 0 (t)/{l + a) 


The dimensionless form of Equation (2.15) and (2.16) are derived in Appendix A. The 
pre-ablation computations are straight-forward once the time at which ablation begins 
is determined from the condition given by (3.9). The solution gives the ablator surface 
temperature ratio, (T, — Tj,)/(T 0 — Tj), the dimensionless thermal thickness, 6 and the outer 
structure (shield) temperature ratio, T w /T*. During ablation the differential equation for 
8 is solved numerically and the shield temperature ratio and dimensionless mass flux axe 
determined as functions of dimensionless time t/t m . 
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The time at which ablation ends ( V = 0) is found by checking the numerical solution 
as it is marched in time. The post-ablation computations are again straightforward and 
give the ablator surface temperature ratio, dimensionless thermal thickness and shield 
temperature ratio as functions of dimensionless time t/t m . 

The characteristic heating rate function F{t/t m ) used as an example in this analysis is 
shown in Figure 3.1. Its shape is representative of a family of curves corresponding to an 
equilibrium glide re-entry trajectory. 

The function F(t/t m ) is computed using the hypersonic approximation to the stagna- 
tion point heating rate for a spherical surface 4 . 




0.763Pr-°- 6 



(3.13) 


where Pr is the Prandtl number, j. i e and p e are the viscosity and density coefficients at the 
boundary layer edge, dUJds is the velocity gradient at the stagnation point, h aw is the 
adiabatic wall enthalpy (equivalent to the total enthalpy, H 0 as used in this paper), and h w 
is the actual wall enthalpy. The velocity gradient is determined from modified Newtonian 
theory^ 5 ’ 6 ! 


ds 



(3.14) 


where r yy is the radius of curvature at the stagnation point, P e is the stagnation point 
pressure, and Poo is the free-stream static pressure. Using the hypersonic approximation, 
h w « h aw , equation (3.13) is simplified and expressed as a function of free-stream density 
and velocity, both functions of time for a known trajectory. Thus, q 0 — q 0 (t) from which 
the heating function F(t/t m ) may be determined. 

Figure 3.2 depicts the total heat input function G(t/tm ) which is determined by nu- 
merically integrating (3.13). Both functions are computed from a trajectory of discrete 
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points and then smoothed by interpolation. 

The characteristic relation between the ablation and radiation parameters is deter- 
mined for this trajectory from equation (3.9) and the condition F{t a /t m ) = R, from 
equation (3.6). This is shown in Figure 3.3. As mentioned above, values of A and R 
which lie above the curve characterize suitable material properties (ablator and shield) 
and insure that T w will remain below the critical limit T* during the pre-ablation period. 
During ablation, the ablator must also satisfy (3.11) to assure this condition for the entire 
trajectory. 

In addition to the ablation and radiation parameter curve generated, several other 
quantities are governed by the value of the ablation parameter. Figure 3.4 is a characteristic 
curve for this example of the relation between the ablation parameter and the time at which 
ablation begins during the trajectory. The curve is a function of the ablation parameter 
and shows that the dimensionless ablation time (t a /t m ) increases with larger values of the 
ablation parameter hence the onset of ablation is delayed as the ablation parameter is 
increased. The values of the dimensionless heat flux and total heat input at the time of 
ablation are also determined by the ablation parameter, although the relation is an indirect 
one since these values are fixed by the ablation time. The characteristic curves are shown 
in Figures 3.5 and 3.6. 

Two other quantities indirectly dependent on the ablation parameters are the the ther- 
mal thickness when ablation begins, plotted as a parameter in dimensionless form in Figure 
3.7, and the parameter ft used to simplify equation (2.15) and (2.16). The parameter ft is 
defined in equation (2.18) and is used to simplify the form of the dimensionless equations 
solved during the period of ablation. The distinguishing feature of the parameter ft is 
that it does not vary very much from a value of one for any given value of the ablation 
parameter (Figure 3.8). This allows further simplification of the governing equations and 
may be used to obtain an analytical solution when a quick zero-order analysis is desired. 
In order to visually demonstrate the concepts discussed in this section, the governing equa- 
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tions derived in Section 2.0 were solved for several ablating materials including water as a 
standard for comparison. The results are presented and discussed in the next section. 

4.0 Results and Discussion 

The one-dimensional thermal analysis presented in the previous two sections led to 
the derivation of two coupled equations to be solved numerically for a given trajectory. 
In this section, the results of such a numerical simulation are presented and discussed. 
The thermal properties of the four ablative materials 7 used as an example axe presented 
in Table 4A. The quantity HC is the threshold enthalpy required for vaporization and 
replaces the quantity cj,(T 0 — 7&) in equation (2.15). It is assumed that the outer shield 
structure is limited to T* = 2000 AT and has an emissivity e = 0.8. The efficiency factor is 
assumed constant at tj = 0.5 for simplicity. Table 4A also includes the value of the ablation 
and radiation parameters for each material as well as the time at which ablation begins 
and the time at which ablation ends. The total mass released by the ablative material per 
unit area of shield to be protected is also tabulated. This is determined by numerically 
solving the coupled dimensional equations (2.15) and (2.16). 

For comparison, the performance of each material may be compared with that of water. 
Water is used as a standard because it has some very favorable properties as an ablator. 
Its relatively high heat capacity allows a high threshold enthalpy to add to a relatively high 
enthalpy of vaporization. Combined with a low vaporization temperature, these properties 
make water an almost ideal ablator for this application. The characteristic enthalpies allow 
it to absorb large quantity of thermal energy without excessive mass loss, while the low 
vaporization temperature insure that a sufficient amount of vapor will be released for the 
transpiration mechanism to achieve maximum effectiveness. 

The ablator surface temperature ratio for the pre-ablation and post-ablation periods is 
shown in Figures 4.1 and 4.2. The curves for lithium flouride lie apart from the clustering 
of the other three primarily because of its much greater characteristic ablation tempera- 
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ture T a . It also has a greater conductivity than the other ablators and thus its ablation 
parameter places it in the undesirable region of the ahlation/radiation parameter curve 
(Fig. 4.6). 

Figure 4.3 is a time history of the dimensionless thermal thickness 0/ 0 a . The results for 
teflon, polyethylene, and lithium flouride reflect monotonically increasing values, while the 
RTV-602 ablator actually depicts a condition in which the recession rate of the ablating 
surface exceeds the rate at which thermal energy is conducted into the material. This 
causes a net decrease in the thermal thickness during ablation (RTV-602 has the lowest 
conductivity). The primary reason behind this is the low energy-absorption capacity per 
kilogram of material. The low vaporization enthalpy results in massive ablation, releasing 
a large amount of gaseous products. Although this effectively protects the external shield 
by transpiration, the total mass loss is excessive. 

The dimensionless mass flux is shown in Figure 4.4. In this case, teflon, RTV-602, 
and to a slightly lesser degree, polyethylene quickly reach a steady ablation state (V=l) 
and remain there until ablation ends. Lithium flouride does not reach a steady state and 
ablates only for a relatively short period of time. The reasoning behind the material’s poor 
performance may be traced to its high vaporization temperature and mediocre energy- 
absorption capacity as it ablates. Despite these results, lithium flouride does maintain 
the outer shield at or below its maximum temperature. Total mass loss is also less than 
the amount of water required, but teflon and polyethelene are slightly lower. As seen in 
Figure 4.6 the material is borderline with regards to the characteristic ablation/radiation 
parameter curve. One significant feature which disqualifies lithium flouride as a good 
choice is that during the post-ablation period, the thermal thickness increases very rapidly 
(Figure 4.3), a condition which would allow and excessive amount of thermal energy to 
penetrate the virgin material. For a finite ablation layer, this would require a prohibitive 
amount of material to protect the interior of the vehicle. 

The shield temperature ratio is plotted in Figure 4.5 for the materials discussed. All 
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four materials performed well in this respect, with RTV-602 keeping the lowest profile. As 
previously indicated, this is due to its low total enthalpy of ablation (HC + HL), which 
allows it to ablate to a higher degree than the other materials, thus demonstrating the 
effectiveness of the transpiration cooling effect. 

Based on the approximate analysis presented, the prime candidate for an ablation 

material in this application should have a high absorptivity coefficient to maximize the 

radiation parameter for a given outer shield. A combination of thermal properties should 

* 

yield an ablation parameter that will place the material as close as possible to the char- 
acteristic R-A curve. This will delay the onset of ablation while still keeping the shield 
temperature below T* prior to ablation. The total enthalpy of ablation (HC+HL) should 
be high enough to allow the material to absorb a maximum amount of thermal energy and 
keep total mass loss at a minimum. Vaporization temperature should be low enough so 
that a sufficient amount of ablation products axe released to take advantage of the tran- 
spiration mechanism. The relative importance of the various thermal properties discussed 
will vary with the trajectory and different materials will be suitable for some and unsuit- 
able for others. The key, nevertheless, lies with the dimensionless parameters A and R 
derived in this analysis as well as constraints placed on the enthalpy of ablation for any 
given trajectory. 

These conclusions are supported by the results of the numerical simulation presented 
above, which clearly favor materials with a combination of thermal properties similar to 
teflon or polyethylene. On the basis of this approximate analysis, these materials per- 
formed according to what was hoped, maintaining the outer shield structure at or below 
its maximum temperature while utilizing a minimum amount of mass per unit area pro- 
tected. 

Although simplified to a one- dimensional analysis, the thermal protection mechanism 
proposed in this paper appears to have promising applications. Besides the hypersonic 
trans-atmospheric vehicles assumed, other aerospace structures would benefit as well, such 
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as inter-planetary capsules upon return to earth, atmospheric orbital transfer vehicles, and 
any other applications where excessive aerothermodynamic heating is of critical concern. 

The next step toward possible application of this concept is to extend the analysis to 
more realistic geometries and to include the coupling effects of the flowfield with the ther- 
modynamics of ablation and conduction. Finite difference analysis to model the unsteady 
heat conduction taking place both in the shield and the ablator material is required to 
more accurately determine how far the thermal energy penetrates into the virgin material. 
In addition, families of curves corresponding to various heating rate histories need to be 
categorized to determine the parameters of most importance. The characteristic relations 
between the various parameters introduced will vary with trajectory, but it is expected 
that the conclusions drawn will be universally valid. 
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5.0 Conclusions 

A new method of thermal protection for a trans-atmospheric vehicle subject to high 
heating rates has been presented. By placing an ablating material behind a fixed-shape 
porous shield, the effectiveness of transpiration cooling may be achieved with the sim- 
plicity of a passive means of thermal protection. The simplified one-dimensional analysis 
treated above applies only in stagnation point regions but may be extended to include 
more complete geometries. Non-dimensional forms of the governing equations has led to 
the derivation of the controlling ablation and radiation parameters. These parameters 
axe characteristic of the thermal properties of the ablation and shield materials for a given 
heating rate history and determine their potential to protect the surface from extreme tem- 
peratures during the pre-ablation period. In addition, the ablating material is required 
to absorb a sufficient amount of thermal energy (characterized by the heat of ablation or 
vaporization) if the outer shield temperature is to be kept below a specified maximum 
during periods of extreme heating when ablation takes place. It was also found that a 
relatively low vaporization temperature is favorable in order to allow the release of gaseous 
products and take maximum advantage of the transpiration mechanism. 

A non-dimensional numerical simulation was performed using four typical ablator ma- 
terials and trajectory corresponding to an equilibrium-glide re-entry into earth’s atmo- 
sphere. The results indicate that ablation materials with thermal properties similar to 
teflon meet the criteria necessary for a successful application of this method based on the 
one-dimensional approximation. The use of this simplified approach indicates that the 
method has application potential; however the simple model presented requires further re- 
finement to more accurately simulate unsteady conduction processes. Ultimately, detailed 
computational results should couple the external flowfield with the internal thermody- 
namics. Experimental verification is also desired prior to prototype implementation of the 
ablat ion/ transpiration method of thermal protection. 


* 
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7.0 Appendix A 


During ablation equations (2.15) and (2.16) of Section 2.2 must be solved to give the 
mass flow rate, dm/dt , and integral thermal thickness, 9 . The equations are 

l~a q °^ ~ \f h ( Ta ~Tb) + L + + pbCb(T a - T b )^ (A.1) 

and 


a r r i a tt 1 d/7z h(T a -T b ) 

+ r+^J it + — r - 1 


(A. 2) 


These are reduced to non-dimensional form by introducing the dimensionless variables 
r, S, V and e defined a is 

time 


r = t/t a - 1 

thermal thickness 


(A. 3) 


mass flux 


6 = 0/0 a 


[ c b(T a — T),) + L + 

(*) r]H 0 }dm/dt 

i 


<Io(t ) 


and the parameter 


(4.4) 


(4.5) 


I 

I 


CbjTg — Tt) 

c b (T a -T b ) + L + rjH 0 
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(4-6) 


i 



Equations (Al) and (A2) become 


and 



(A.7) 


V(1 ~ «) + 


gg 

go<$ 


1 


where = d6/dr and the relations 


(A.8) 


Q a = p b c b (T a - T b )6 a 


g 0 = k b (T a — T b )/9 a 


and 

0a = Qa^a/Qa 

have been used. Eequations (A.7) and (A. 8) are coupled non-linear differential equations 
for 6 and V with all other quantities including q 0 (r), Q a (r), and e(r) known. A simple 
approximate solution for 6 and V can be found for small values of e as follows: 
Elimination of V from equations (A.7) and (A. 8) gives 



the solution for e = 0 is 


(A.9) 


S 1 = 1 + 2 p a T 


(A.10) 
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Similarly, dimensionless mass flux V is given for e = 0 by 


V = 1 ^ — r (-4.11) 

(1 + 2 0*t) 2 

These may be extended to first order in e by returning to equations (A. 7) and (A. 8). 
Thus, to 0(e), the solution is 

8 2 = 1 + 23 a r - I* 2f3 a [(q 0 (Q/q a )(l + 2&0* - l]e(C)d( (A.12) 

Jo 

and 

y _ i _ (gq/g») , L _ (gq/go) 

( 1+2 P a )T? [ (1 + 2 Pat)~K 

tA f [«.(()/«.)(! + 2AO* - 1]«(CK (j4.13) 

(1 + 2p 0 r )2 jo 

The results presented in Section 4.0 of this paper were obtained by direct numerical 
integration of equations (A. 7) and (A.8). The zero-order approximation was used to derive 
the relation given by the inequality (2.17) or (3.12) which imposes a condition to be satisfied 
by the ablating material in order to avoid excessive temperatures at the outer surface. 
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Fig. 1.1 Stagnation Point Geometry. 




qo(t)/qm 






Fig. 3.3 Characteristic Curve: Radiation and Ablation Parameters. 



Fig. 3.4 Ablation Time vs Ablation Parameter 
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Table 4A : Thermal Properties of Four Ablators and Water 


Material: 

Teflon 

Polyeth. RTV-602 Lith. 

Flrd. 

H 2 O 

Absorptivity 

0.667 

0.667 

0.667 

0.667 


Ablat. Temp. T a (K) 

813.0 

608.0 

707.0 



Back Temp. T b (K) 

298.0 

298.0 

298.0 

298.0 

298.0 

Density (kg/m^) 


950.0 

1000.0 

2640.0 

1000.0 

Heat Cap. (J/kg/K) 

1108.5 

3516.3 

2593.82 

2025.8 

4200.0 

Conduct. (W/m/K) 


2.320 

0.278 

5.749 

0.600 

HC (threshold: J/kg) 

4.93xl0 4 

8.59x1 0 4 

8.54xl0 4 

3.83x10® 

3.15x10® 

HL (ablation; J/kg) 

6.56x10® 

4.65x10® 

4.06x1 0 4 

7.74x10® 

2.26x10® 

ABL. Parameter.A 

0.000282 

0.00106 

0.000172 

0.122 

2.03x10-® 

RAD. Parameter, R 

0.5143 

0.5143 

0.5143 

0.5143 

0.5143 

t a (sec) 

84.9 

127.5 

72.9 

460.5 

37.11 

t e (sec) * 

2279.2 

2259.6 

2261.0 

1892.0 

2282.4 

Mass Req'd (kg/m^) 

120.9 

146.4 

418.7 

169.0 

210.1 


* Total time for heat pulse is approximately 2300 seconds. 


30 









(Ts-Tb) 

(Ta-Tb) 



Fig. 4.1 Ablator Surface Temperature Ratio (Pre- Ablation). 


(Ts-Tb) 

(Ta-Tb) 



Fig. 4.2 Ablator Surface Temperature Ratio (Post- Ablation). 
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Radiation Parameter, 



06 



Ablation Parameter, \[A~. 


Fig. 4.6 Characteristic Parameter Curve with Four Ablators. 
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